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2.  Nonlinear  Modulational  Stability  of  an  Electromagnetic  Pulse  in  a 
Two-Component,  Neutral  Plasma,  by  D.J.  Kaup  and  Ronald  E.  Kates 
[Accepted  by  j.  Plasma  Physics] 

This  manuscript  describes  how  to  correctly  calculate  the 
nonlinear  coefficients  in  the  case  of  an  electromagnetic  pulse 
propagating  in  a  two-component  plasm, a.  We  also  demonstrate  that 
other  values  given  in  the  literature  are  incorrect.  We  correct 
the  predictions  for  such  electromagnetic  propagation  and  discuss 
the  astrophysics!  consequences. 

The  above  results  have  been  presented  as  short  talks  at  two 
meetings,  the  APS  pi  asm a  physics  meeting  in  Nov.  1988  and  at  the 
Grossman  general  relativity  meeting  in  Australia  in  Aug.  1988. 


3.  The  Influemce  of  an  Amibient  Magnetic  Field  on  the  Nonlinear 
Modulational  Stability  of  Circularly  Polarized  El ec tromagne t i c 
Pulses  in  a  Two-Component,  Neutral  Plasma,  by  Rona’d  Kates  and  D.J. 
Kaup  [accepted  by  J.  PI  asms  Physics]. 

Here  we  demonstrate  that  in  the  1  i  mi  it  of  a  strong  magnetic 
field,  the  modu 1  at i onal 1 y  stable  case  in  #2  above  will  become 
modu 1  a t i on  a  1 1 y  u  n  s  t  ab 1 e . 

The  following  work  is  in  various  stages  of  being  completed. 
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4.  The  Modulational  Instability  of  a  Relativistic  Electromagnetic 
Pulse  in  a  Two-Componen t ,  Neutral  Plasma,  by  Ronald  E.  Kates  and 
D.J.  Kaup  (in  preparation). 

This  problem  has  finally  simplified.  In  particular,  the 
longitudinal  electric  field  is  found  to  be  a  dominate  component 
in  this  instability.  And  we  can  now  extend  the  results  on 
modulational  instability  from  the  nonre  1  at  i  v  i  st  i  c  domiain  up  into 
the  relativistic  doma i n . 

5.  A  Singular  Perturbation  Expansion  for  the  Cavity  Modes  of  a 
Micron-Size  Spherical  Droplet,  by  D.J.  Kaup  [in  preparation] 

This  paper  presents  a  rapid  proceedure  for  generating  analytic 
solutions  for  cavity  modes  of  such  droplets,  although  we  have 
only  calculated  the  e i genf r equenc i es  so  far,  the  potential  is 
now  here  to  an  a  1 y  t i c  a  1  1 y  evaluate  non  1  inear  interactions  in  such 
droplets,  and  under  a  verity  of  conditions. 
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We  study  nonlinear  effects  including  possible  modulational  instability  of  an  intense 
electromagnetic  pulse  propagating  through  a  fully-ionized,  unmagnetized  plasma..  (The 
pulse  is  assumed  to  be  strong  enough  to  accelerate  particles  to  weakly,  but  not  fully, 
relativistic  velocities.)  ^The  envelope  is  shown  to  evolve  over  long  time  scales  in  general 
according  to  a  vector  form  of  the  well-known  cubic  nonlinear  Schroedinger  (NLS)  equation. 
Three  distinct  nonlinear  effects  contribute  terms  cubic  in  the  amplitude  and  thus  can  be  of 
comparable  magnitude:  ponderomotive  forces,  relativistic  corrections,  and  harmonic  gen¬ 
eration.  In  contrast  to  previous  work,  our  calculation  takes  all  three  effects  into  account. 
Integrability  and  modulational  stability  of  any  given  system  are  shown  to  depend  on  polar¬ 
ization,  frequency,  composition,  and  temperature,  and  these  dependences  are  given.  Finite 
temperature  effects  are  considered  for  completeness;  they  are  qualitatively  important  only 
near  a  critical  frequency  (<  3/2  u>p)  just  above  the  plasma  frequency  for  which  the  group 
velocity  is  comparable  to  either  (or  both)  of  the  sound  velocities.  In  the  special  case  of  a 
cold  positron-electron  plasma  -  contrary  to  the  predictions  of  Chian  and  Kennel  [1  ]  -  the 
model  is  strictly  mod ulationally  stable  for  both  linear  and  circular  polarization;  The  results 
have  important  implications  for  pulsar  micropulse  observations  and  possible  technological 
applications. 
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The  propagation  of  a  relativistically-strong,  modulated  electromagnetic  pulse  in  a 
plasma  is  more  complicated  than  appears  at  first  glance.  In  the  case  of  uniform  propagation 
(i.e.,  no  modulation),  some  exact  solutions  (for  circular  polarization)  are  indeed  known  [2] 
which  are  described  by  a  nonlinear  dispersion  relation.  Chian  and  Kennel  [1]  pointed 
out  the  possible  application  to  pulsar  micropulses  and  obtained  a  scalar  NLS  equation 
by  applying  the  method  of  Karpman  and  Krushkal  [6]  (which  is  based  on  the  Whitham 
averaging  method  [7]  )  to  this  dispersion  relation.  The  problem  was  again  considered 
by  Mofiz  and  collaborators  [3-5]  Unfortunately,  the  resulting  NLS  coefficients  in  these 
papers  are  incorrect,  and  the  force  equations  in  [1]  are  violated  at  leading  order.  (The 
necessity  of  paying  special  attention  to  the  longitudinal  force  equations  was  noticed  in  [3].) 
The  problem  evidently  contains  subtleties  and  requires  the  careful  application  of  singular 
perturbation  methods. 

The  scalar  NLS  equation  is  given  in  standard  dimensionless  form  by 


*  ±2$*  q 2 


(1) 


The  upper  sign  corresponds  to  modulational  instability.  As  we  shall  see,  the  propagation 
of  the  envelope  of  an  electromagnetic  pulse  in  an  unmagnetized,  neutral,  two-component 
plasma  in  the  weakly  nonlinear  case  is  governed  in  general,  not  by  (1),  but  by  a  vector 
form  of  the  NLS  equation.  By  weakly  nonlinear  we  mean  that  the  transverse  electric  field 
accelerates  one  or  (possibly)  both  components  to  weakly  relativistic  velocities  within  one 
cycle,  that  is, 
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In  the  general  case,  two  NLS  amplitudes  will  be  present,  corresponding  to  the  two  po¬ 
larization  states.  However,  these  polarization  states  will  be  coupled  by  the  nonlinearity. 
Consider  a  system  of  the  form 

idrQi  +  Dd%q\  +  q\{C\q\q *  +  C29292)  ~  0  (3a) 

idr<j2  +  Dd^qi  +  92(^192^2  ■+■  C2919*)  =  0  (3 b) 

In  the  special  case  where  C\  and  C2  are  equal  (or  can  be  transformed  into  this  form),  the 
equations  are  known  to  be  integrable  [8].  However,  in  general,  the  vector  NLS  equation 
will  be  nonintegrable. 

The  purpose  of  this  paper  is  to  present  an  accurate  evaluation  of  the  leading  non¬ 
linear  coefficients,  because  they  are  crucial  in  determining  the  qualitative  behavior  of  the 
envelope,  in  particular  whether  or  not  modulational  stability  and  soliton  type  solutions 
will  occur.  The  reasons  for  performing  this  calculation  are  several  fold.  For  example,  with 
the  advent  of  current  technological  capabilities,  one  could  be  interested  in  how  nonlinear 
effects  will  couple  different  polarizations,  particularly  if  one  intends  to  have  each  polar¬ 
ization  mode  carry  different  bits  of  information.  Also,  intense  signals  in  a  rarified  plasma 
could  have  critical  information  scrambled  by  the  nonlinear  effects.  It  may  even  be  possible 
to  test  our  theory  in  astrophysical  systems,  in  particular  by  observing  pulsar  micropulses. 
This  point  will  be  discussed  by  us  in  a  separate  paper  [9]. 
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In  any  singular  perturbation  method,  one  seeks  uniformly  valid  asymptotic  expansions 
[10]  with  respect  to  a  small  parameter,  which  we  will  call  e  (in  this  problem  e  is  related 
to  the  typical  amplitude  of  the  vector  potential).  Here,  the  domain  of  uniform  validity  is 
of  order  1/e  cycles  in  space  and  1/e2  cycles  in  time.  The  singular  perturbation  method 
which  we  will  apply,  known  as  the  method  of  multiple  time  scales,  has  been  described  in 
several  textbooks  (see  for  example  Van  Dyke  [11]  or  Nayfeh  [12])  and  has  been  applied  to  a 
wide  variety  of  problems.  Since  we  are  interested  in  the  evolution  of  a  pulse  with  arbitrary 
modulation,  it  will  not  be  possible  to  assume  dependence  on  only  a  single  independent 
variable  (although  it  is  possible  to  verify  consistency  of  our  results  for  particular  cases 
with  the  work  of  Kozlov,  et.  al.  [13],  in  which  fully  nonlinear  solutions  depending  on 
one  independent  variable  were  found).  Instead,  we  allow  the  functions  in  this  problem  to 
depend,  not  only  on  a  “fast”  (phase)  variable,  but  also  on  “slow”  time  and  space  variables. 
These  will  be  chosen  to  make  the  effects  of  nonlinearity  enter  the  calculation  at  the  same 
order  as  the  effects  of  linear  spreading  of  the  wave  packet. 

Consider  an  envelope  of  length  L  of  waves  with  wavelength  near  A.  If  c  =  ^“l^l  ~  r, 
then  nonlinear  effects  will  typically  act  on  a  timescale  To  =  A/(e2c)  comparable  to  the 
timescale  for  the  spreading  of  packet  due  to  linear  dispersion.  By  routine  dimensional 
analysis  it  is  possible  to  see  that  ponderomotive  forces,  relativistic  corrections,  and  har¬ 
monic  generation  might  have  comparable  effects  on  this  timescale.  Our  calculation  shows 
that,  in  general,  all  three  effects  indeed  play  a  role. 

The  ponderomotive  force  depends  quadratically  on  the  amplitude  and  leads  to  a  slowly 
varying  longitudinal  field,  corresponding  physically  to  radiation  pressure,  which  leads  to 
slow  longitudinal  motions  and  modifies  the  background  density.  At  cubic  order,  the  modi¬ 
fication  of  the  background  couples  back  to  the  fundamental.  This  contribution  dominates 
at  the  longest  wavelengths  (c7k2  u>2).  (In  fact,  for  finite  temperature,  there  is  a  sin¬ 

gularity  at  a  critical  frequency  just  above  the  plasma  frequency.)  At  shorter  wavelengths 
(c2*2  ;>  a>2),  the  other  two  contributions  dominate  in  an  electron-positive  ion  plasma, 
while  all  three  effects  are  important  in  an  electron-positron  plasma.  Independently  of  the 
composition,  temperature,  or  frequency,  the  effects  of  harmonic  generation  are  identically 
zero  for  circular  polarization  -  but  in  general  only  for  circular  polarization. 

Chian  and  Kennel  [1]  first  proposed  that  modulational  instabilities  in  an  electron- 
positron  plasma  would  reshape  pulsar  micropulse  signals.  Unfortunately,  they  omitted 
two  ’’parametric”  sources  of  nonlinearity  (harmonic  generation  and  ponderomotive  effects) 
and  also  obtained  an  incorrect  value  for  the  linear  dispersion  coefficient  of  the  NLS.  As  we 
shall  see,  the  special  case  of  an  electron-positron  plasma  exhibits  several  peculiarities,  such 
as  accidental  (near)  cancellations  in  certain  limits.  One  can  obtain  the  correct  nonlinear 
coefficient  only  by  considering  all  three  sources  of  nonlinearity  right  from  the  beginning. 
When  all  terms  are  included,  we  find  that,  at  least  in  the  nonrelativistic  limit  (their  Eq. 
(10)),  their  claims  concerning  modulational  instability  are  reversed  -  both  for  linear  and 
for  circular  polarization.  A  cold  positron-electron  plasma  is  modulationally  stable  for  both 
circular  and  linear  polarization.  However,  if  ions  are  present,  modulational  instability  is 
possible,  as  we  shall  see.  Implications  for  pulsar  signals  will  be  briefly  discussed  below. 

Our  model  comprises  a  neutral,  fully-ionized  plasma,  consisting  of  singularly-ionized 
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atoms  (of  charge  +e  and  mass  M)  and  electrons  (of  charge  — e  and  mass  m);  the  plasma 
is  modeled  as  a  fluid  with  constant  but  distinct  electron  and  ion  temperatures  (possibly 
zero).  The  fluid  is  characterized  by  the  electron  density  n,  the  electron  velocity  v ,  the  ion 
density  N  and  the  ion  velocity  V.  (Although  we  will  refer  to  the  positive  component  as 
”the  ions,”  our  model  of  course  also  describes  a  positron-electron  plasma  if  m  =  M.)  The 
plasma  is  further  characterized  by  a  vector  potential  A  and  an  electrostatic  potential  <f>. 
The  functions  n ,  v  N.  V ,  A,  and  4>  are  assumed  to  depend  on  t  and  z,  but  not  on  x  and  y. 
It  is  convenient  to  write  the  continuity  and  Euler  equations  for  the  electrons  in  the  form 


dtn  +  V  •  (ntT)  =  0 

d(myv  —  -A)  —  eVd  +  -E^-jV/VA*  —  mc^V  Inn  =  0 
c  c 

where  cs  is  the  electron  speed  of  sound,  m  is  the  electron  mass,  and 


(4a) 

(46) 


d=dt  +  v-$  (4c) 

7  =  (l-vVc2)-1/2  (4u) 

For  the  ions,  the  respective  equations  are 

dtN  +  V  •  (NV)  =  0  (5a) 

D(MVV  +  -A)  +  eV<t>  -  -ELi V,VAe  -  MC ?  •  V  In  JV  =  0  (56) 

c  c 

where  C ,  is  the  ion  speed  of  sound,  M  is  the  ion  mass,  and 

D  =  dt  +  V  •  V  (5c) 

r  =  (i  -  v2/c2)-1/2  (5  d) 


Here,  the  electron  speed  of  sound  c ,  is  assumed  constant  and  is  related  to  the  temperature 
Te  by  c]  =  'yeTe/m>  and  similarly  for  the  ions. 

Of  course  Maxwell’s  equations  must  hold: 

VxB  =  ~(nv~NV)+  -dtE  (6a) 

c  c 

V  -E  =  4ire(N  -  n )  (66) 


with 

B  =  V  x  A 

E  =  -Vd  -  -dtA 
c 
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(6c) 

(6*0 


and  where  we  impose  the  radiation  gauge 

V-A  =  0  .  (6e) 

This  is  a  nonlinear  system  of  equations  is  to  be  solved  for  the  evolution^of  some  initial 
data  given  on  the  interval  — oo  <  z  <  oo  for  the  functions  (n,  t7,  N ,  V,  A ,  <£),  satisfying 
the  usual  initial-value  constraints.  (The  constraints  will  be  satisfied  automatically  at  later 
times).  We  are  interested  in  a  particular  class  of  solutions,  namely  those  which  correspond 
to  the  physical  idea  of  a  slowly  modulated,  weak,  nearly  sinusoidal  disturbance  about  a 
uniform  quiescent  medium  (which  is  an  obvious  exact  solution  of  Eqs  (4-6)).  We  assume 
that  well-posed  initial  data  can  be  given  which  evolve  in  this  way. 

To  make  these  assumptions  precise,  consider  first  the  electron  density  and  the  vector 
potential,  which  are  assumed  to  have  asymptotic  expansions  of  the  form 

n  =  no  +  eri]  +  €* ri2  +  e3n3  +  . . .  (7a) 

A  =  tA\  +  t*  A2  +  c3A3  +  . . .  (76) 

with  respect  to  the  positive  dimensionless  parameter  e  1,  where  the  first  term  n0  in  (7a) 
is  a  constant.  Since  the  velocities  will  also  scale  with  e,  it  is  evident  that  for  sufficiently 
small  e  the  motion  is  only  weakly  relativistic  and  thus  we  will  be  able  to  expand  the 
relativistic  factors  7  and  T  about  unity. 

Similarly,  the  other  variables  are  expanded  as  follows: 


N  ■=  n0  +  eN\  +  N 2  A  f3 N3  4-  . . . 

(7c) 

V  =  €V~i  +  e2V2  +  t3V3  +  •  • . 

(7  <0 

V=  tVi  +e2F2  +  e3F3  +  ... 

(7a) 

<t>  =  t<f>  1  -1-  e2  $2  -f  e3  <f>3  +  •  •  • 

(7/) 

Note  that  no  =  const,  No  =  no,  and  all  other  zeroth-order  quantities  vanish.  (Except 
where  otherwise  stated,  numerical  subscripts  indicate  the  formal  order  in  e.) 

Now,  it  is  routine  to  work  out  the  linearized  solutions  of  (4)-(6)  corresponding  to  a 
sum  over  Fourier  components  of  transverse  electromagnetic  EM  waves.  (By  ’’linearized” 
solutions,  we  mean  solutions  obtained  by  simply  truncating  (4)-(6)  at  order  e.)  A  typical 
Fourier  component  of  the  vector  potential  perturbation  behaves  in  this  linearized  theory 
like  a(u>)el(kz~ut\  where  k  and  u  are  related  by  a  dispersion  relation  to  be  given  below. 
Here  we  are  interested  in  the  nonlinear  theory  in  which  different  Fourier  components  are 
coupled.  In  order  to  study  the  evolution  of  a  slowly  modulated  beam,  we  assume  that  the 
first-order  perturbation  of  the  vector  potential  takes  the  form 

A,  =  aei{kx~wt)  +  c.c.  (8) 

where  “c.c.”  means  the  complex  conjugate,  k  is  the  wave  number,  u>  is  the  frequency.  The 
functions  ax  and  ay  (where  the  subscripts  here  refer  to  the  component)  are  not  strictly 
constant,  but  depend  on  the  slow  variables 

Z  =  ez,  T  ~  et,  r  =  e2<  (9) 
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With  the  introduction  of  slow  coordinates,  we  employ  the  following  standard  procedure: 
We  collect  terms  in  all  equations  according  to  their  formal  order  of  magnitude.  (“Order” 
always  means  with  respect  to  e.)  We  then  solve  order  by  order.  Every  derivative  of  a  slowly 
varying  function  is  down  by  one  or  two  orders.  For  this  reason,  the  transverse  equations 
at  leading  order  are  the  same  as  those  of  the  truncated  linearized  theory.  (Not  so  for  the 
longitudinal  equations,  as  we  shall  shortly  see.)  Slow  derivatives  of  the  lower-order  terms 
will  appear  at  higher  orders.  According  to  standard  singular-perturbation  procedure  [12], 
we  seek  solutions  of  the  higher-order  equations  which  satisfy  the  requirement  of  uniform 
validity,  not  only  over  a  cycle,  but  also  over  the  slow  scales.  (In  general,  not  all  of  the 
higher-order  equations  need  to  be  solved  explicitly,  but  instead  provide  integrability  con¬ 
ditions  which  affect  the  lower  orders.)  In  this  way,  the  evolution  of  the  slowly  varying 
functions  will  be  determined. 

With  this  strategy  in  mind,  we  now  solve  the  first-order  equations.  The  longitudinal 
components  are  identically  zero  or  can  be  chosen  to  vanish  without  loss  of  generality.  The 
transverse  components  can  be  treated  as  a  coupled  linear  homogeneous  system  of  the  form 

m i)  =  o  .  (io) 


At  this  order,  we  choose  the  EM  mode  to  dominate  and  take  the  solution  to  be  of  the  form 


t>i  = 


^>1  =  0  =  na  —  N\ 


Au  =  vi. 


Vi, 


(Ha) 

(116) 

(11c) 


Note  that  the  velocities  and  the  vector  potential  are  purely  transverse  in  first  order.  The 
first-order  part  of  (6b)  guarantees  quasineutrality  at  first  order.  The  z-component  of  (6a) 
is  then  satisfied  at  first  order  automatically. 


The  dispersion  relation  is 

w2=wj  +  c2fc2  (12) 

where 

(i  +  jf) 

The  slow  dependence  of  a  is  of  course  as  yet  undetermined 

As  we  have  seen,  the  longitudinal  components  of  the  velocities  vanish  in  first  order. 
However,  for  a  system  with  slowly  modulated  waves,  the  third-order  part  of  the  longitudinal 
component  of  the  force  equations  (4b)  and  (5b)  leads  to  a  nonvanishing,  slowly  varying 
longitudinal  component  of  the  velocities  at  second  order.  This  happens  even  in  a  positron- 
electron  plasma,  despite  the  fact  that  the  third-order  longitudinal  electric  field  (second- 
order  electrostatic  potential)  vanishes.  Note  that  in  Ref.  [1],  the  longitudinal  component 
of  the  force  equations  (their  Eq.  (5))  is  violated  in  leading  order. 


(13) 
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At  second  order  there  are  three  possible  types  of  terms  with  different  phase  depen¬ 
dence,  which  we  will  discuss  in  turn.  Decomposing  the  second-order  electron  density  n 2 

n2  =  n<0)  +  [r4,)ei(fc*-wl)  +  c.c.]  +  [n^e2**—"0  +  c.c.]  ,  (14) 

we  note  that  the  first  term  n2(0)  corresponds  to  a  second-order  correction  to  the  back¬ 
ground  density.  A  similar  decomposition  will  be  made  for  the  other  variables.  From  here 
on  superscripts  (0),  (1),  (2),  etc.  refer  to  coefficients  of  terms  with  corresponding  phase 
dependence.  Terms  such  as  be  refered  to  as  “DC”  terms.  Terms  such  as 

*  f  2) 

correspond  to  a  correction  to  the  fundamental.  Terms  such  as  correspond  to  the  har¬ 
monic,  which  in  general  arises  due  to  nonlinearity,  although  in  specific  cases  it  may  vanish 
identically. 

At  all  orders  p  >  1,  the  fundamental  mode  is  described  by  an  inhomogeneous  system 
of  the  general  form 

L(4jp)  =  NP(rJ>i,...ipp-i)  -  (15) 

where  L  is  the  linear  operator  of  Eq.  (10),  and  where  Np  is  a  (possibly)  nonlinear  functional 
of  the  lower-order  terms.  If  one  defines  the  scalar  product 


Jo 


de  fg 


of  two  functions  f  and  g  over  a  cycle,  then  Np  must  be  orthogonal  (with  respect  to  (16))  to 
the  general  homogeneous  solution  of  (10).  Otherwise,  secular  terms,  e.g.,  terms  of  the 
form  6e'8  (9  =  kx  —  would  arise.  Secular  terms  would  obviously  violate  the  uniformity 
condition  stated  above,  because  the  pth  order  would  eventually  become  comparable  in  size 
to  the  terms  of  order  p  —  1.  At  second  order,  this  requirement  implies  simply 


(dr  +  Vgdz)  a  =  0 


where 


Note  that  Eq.  (17)  is  the  same  as  the  group  velocity  condition  for  a  wave  packet  in 
the  truncated  linear  theory.  (Nonlinearites  will  play  a  role  in  the  fundamental  at  third 
order.)  Now  that  the  secular  terms  in  the  fundamental  mode  have  been  removed,  then  the 
solution  for  the  secopd-order  fundamental  can  be  taken  to  vanish  without  loss  of  generality 
(a  nonvanishing  contribution  could  always  be  absorbed  into  first  order  )  Thus 

r4,}  =  N2(1)  =  =  0 

rn  (19) 

t£(1)  =  v2(1)  =  A2  =  0  . 


Now  for  the  second  order  DC  terms.  These  are  obtained  by  solving  ^-components  of 
the  Euler  equations  (4b)  and  (5b)  in  third  order  together  with  the  continuity  equations 
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(4a)  and  (5a)  at  second  order. 


n 


(0)  =  jv<«»  = 


e2n0 


Mm  c2v2 


-40)  -40) 


=  z  vgn2 


(o) 


/no 


40)  =  0 


(20a) 

(206) 

(20c) 

(20  d) 


Note  that  in  all  cases  the  plasma  acquires  a  longitudinal  component  of  the  velocity  in 
the  ^-direction,  even  though  <f>^  vanishes  for  an  electron-positron  plasma.  Except  for  the 
electron-positron  case,  the  DC  electric  field  will  induce  charge  separation  over  the  slow 
space  scales  at  <9(e4). 

The  nonlinear  terms  in  the  Euler  equations  (4b)  and  (5b)  also  induce  harmonic  terms 
in  the  solution: 


(2)  _  n0e2 fc2(q  •  a)  ( 4u^  u2\ 

nz  2m  u>2c2( 4u.’2  —  u>2)  y  m  M  J 

*r(2)  =  n0e2k2(a-a)  ( 4u^_  _ 

2  2 M  u>2c2( 4u:2  -  u2)  \M  m  J 


vi2)  = 


wn 


(2) 


kriQ 


V™  = 


wA 


(2) 


4”  =  ^ 


kn0 

2(“~ 
v  m 


)  (a  •  a) 


2c2(4u>2  —  u;2) 


42)  =  0 


(21a) 


1216) 

(21c) 

(21d) 


(2 1  e ) 
(21/) 


Evidently  from  Eqs.  (21),  circular  polarization  (vanishing  a  ■  a)  is  a  special  case  (for  any 
mass  ratm),  because  harmonic  generation  is  absent.  In  this  case,  charge  quasi-neutrality  is 
maintained  at  second  order.  However,  in  general  the  second-order  density  perturbations  of 
electrons  and  ions  are  not  equal:  charge  separation  over  the  fast  scales  is  already  present 
at  second  order! 

The  electron-positron  case  is  again  special  because  the  harmonic,  second  °rder  poten¬ 
tial  perturbation  is  zero.  Otherwise,  a  second-order  electric  field  varying  as  the  harmonic 
will  be  present  in  general  (except  for  circular  polarization). 

We  now  proceed  to  third  order.  It  is  now  possible  to  derive  the  evolution  of  the  slow 
functions  (a.  etc.)  as  a  consequence  of  the  requirement  that  the  third-order  equations 
contain  no  secular  terms,  as  explained  above.  Since  DC  and  harmonic  terms  are  of  course 
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orthogonal  to  the  operator  L  of  Eq.  (15),  it  sufficies  to  consider  the  fundamental.  As  at 
second  order,  one  can  absorb  the  perturbations  in  the  fundamental  mode  of  the  third-order 
electron  density  »4*\  i°n  density  and  potentials  <f>^\  -A^  into  the  first-order  terms 

without  loss  of  generality.  Thus, 

=  ^3(1)  =  =  0,  41)  =  0  •  (22a) 


However,  the:  _  remains  a  third-order  perturbation  of  the  velocities  of  the  form  which  comes 
from  the  relativistic  corrections 


=  2mh*  I2S(“"  '  “»  (22!>) 

^  =  S‘|2i,7'i)  +  !'(i'i)l  (22C) 

The  secular  condition  then  becomes 


2iu.'dra  -f  (c2<9 1  —  dj-)a  —  /no) 


2  -* 

- 


,<2) 


n;  '  Ar, 


+ 


„2,  -2 


mn0 
l 


+ 


(2) 

2 


Mno 


JL  _L  J_ 
m  ‘  M 


+ 


AP 


2c“  £  +  -n- 

TTl  »*] 


[2a(a*  •  a)  -f  a*(a  ■  a)]  =  0 


(23) 


The  origin  of  the  various  terms  in  (23)  is  evident.  The  first  two  terms  would  be  present 
even  in  the  linearized  system  obtained  by  simple  truncation  of  Eqs.  (4-6);  they  represent 
the  effects  of  linear  dispersion.  The  term  containing  arises  because  the  ponderomotive 
force  induces  perturbations  of  the  “background”  density  on  the  slow  scales,  as  discussed 
above.  The  terms  containing  n2  ;  and  N?  '  arise  from  harmonic  generation.  Finally,  the 
last  term  comes  from  relativistic  momentum  corrections  contained  in  the  factors  of  7  and 
T  of  Eqs.  (4)  and  (5). 

Taking  into  account  Eqs.  (17),  (18),  (20),  and  (21),  we  can  rewrite  (23)  in  the  form 


2 iu.'dTa  +  -~-d\a  +  S' (a  ■  3){Ch  +Cr) 

+  ( Cp  4-  2Cfl)a(a*  ■  a)  =  0 


(24) 


In  the  above,  the  subscripts  H  (haimonic),  P  (ponderomotive),  and  R  (relativistic)  indicate 
the  origin  of  the  corresponding  nonlinear  coefficient.  These  nonlinear  coefficients  are  given 
for  zero  temperature  by 


Ch 

CP 


2u’ic2(4oj2  —  u>2) 


e2^p 


4u,’ 


1  _L  1 
rrv?  +  — 


2  m 


At  5 


-L  +  J_ 
m  r  At 


mM 


Mm  c2v2 


(25a) 

(256) 
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(25  c) 


Cr  = 


.2,  ,21,1 

e  up  7^+W 


2c4  J-  +  i 

m  '  M 


If  the  temperature  is  finite,  the  relativistic  coefficient  remains  unchanged,  whereas  Eqs. 
(25a)  and  (25b)  generalize  to 


e2k2u 2 


f  mM  \  f 
\M  +  m  )  1 


+ 


m3(u>2  -  k2c2)  M*(u>2  -  k2C 2) 


+ 


L J 

P 

4£>2 


/  mM  \ 


1 


m2(u>2  —  c2fc2)  M2(u>2-C2k2) 


with 


and 


D 

(25d) 


X>2  =  1  —  T« 


mM 


1 


Cp  =  u> 


4  p  M  +  m  [m(w2 
M  +  m' 


+ 


c2h2)  M(w2-C2/b2)J 


Mm 


(25  e) 


-  C2)  +  m(u2  -  c2) 

(Note  that  the  factor  D2  varies  between  3/4  and  1.)  In  the  high-’wavenumber  limit,  all 
three  coefficients  approach  constants,  whether  or  not  the  temperature  is  finite.  In  the  case 
of  zero  temperature,  as  k  approaches  zero,  the  ponderomotive  term  dominates,  because 
the  coefficient  Cp  diverges  like  k‘K~2\  whereas  the  relativistic  coefficient  Cr  approaches  a 
constant  and  the  harmonic  coefficient  Cr  vanishes  like  k2 .  If  the  temperature  is  finite, 
the  ponderomotive  coefficient  Cp  is  singular  just  above  the  plasma  frequency.  (Since  the 
speed  of  sound  can  never  exceed  l/\/3 c,  this  critical  frequency  lies  below  3/2  u>p). 

In  the  cases  of  linear  and  circular  polarization,  the  vector  NLS  reduces  to  a  scalar 
NLS.  For  circular  polarization,  we  obtain 


2  2 

2  iudra  +  ~^-d2z  a  +  2  (CP  +  2CR)aa2  =  0 


(26) 


where  we  have  set  a  =  a  ■  (x  +  ly).  For  linear  polarization,  we  obtain 


2  iu>drQ-  +  --^-dza  +  (C  r  +  3  Cr  +  Cp)o*  a.2  —  0  (27) 

where  now  we  have  taken  a  =  a  •  z.  Thus,  modulational  instability  for  these  two  cases 
occurs  if  the  combinations  (Cp  -fi  2Cp),  (Cr  +  3 Cr  -f  Cp),  respectively  are  positive  [8]. 

We  now  consider  some  special  cases: 

ELECTRON-POSITIVE  ION  PLASMA:  If  the  plasma  is  composed  of  positive  ions  and 
electrons,  then  m/M  <C  1  and  we  find  (for  the  zero-temperature  case) 


Ch  +  Cr- 


3eX 


2m2c4(4o>2  —  u>2) 


>  0 


CP  +  2Cfi= 


2  2 
e  u>* 

c*m 2 


1  - 


mur 
Mc2k 2 


(28a) 

(286) 
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Except  for  frequencies  very  near  the  plasma  frequency,  both  polarizations  are  evidently 
modulationally  unstable.  The  corresponding  results  for  finite  temperature  are: 


Ch  +  Cr= 


3«X(i  - 

2m2c4(4o;2  —  u>2  —  4  c2h2) 


e2u>2 

Cp  +  2Cr=  — — |- 


0- 


m 


u; 


Afw2(l-C2/c2) 


r*) 


(28c) 
(28  d) 


As  in  the  zero-temperature  case,  both  polarizations  are  modulationally  unstable  -  ex¬ 
cept  ior  a  small  frequency  range  just  above  the  plasma  frequency.  At  finite  temperature, 
the  singularity  in  the  ponderomotive  coefficient  near  the  plasma  frequency  leads  to  corre¬ 
sponding  singularities  for  both  linear  and  circular  polarization;  sufficiently  near  the  plasma 
frequency,  the  system  is  modulationally  stable.  On  the  other  hand,  finite  temperature 
makes  virtually  no  difference  in  the  behavior  of  either  polarization  at  moderate  to  high 
frequencies  -  even  though  the  combination  Ch  -f  Cr  can  become  negative.  In  all  cases, 
the  contribution  of  the  ponderomotive  coefficient  is  down  by  0(m/M )  compared  to  the 
others,  except  near  the  lasma  frequency.  The  positive  (destabilizing)  contribution  of  the 
relativistic  coefficient  outweighs  the  negative  (stabilizing)  contribution  of  the  harmonic 
coefficient. 

ELECTRON-POSITRON  PLASMA:  Next,  let  us  consider  an  electron-positron  plas¬ 
ma,  as  in  Chian  and  Kennel  [1].  Setting  M  =  m  in  Eqs.  (25),  we  find  in  the  zero- 
temperature  case 


Ch+C*=  2 mW  >0 

(29n) 

— e2u>4 

(29b) 

c'  +  2C*=mw<02  2 

Ch  +  ZCr  +  Cp  —  +“,) 

(29c  j 

Since  Ch  +  3 Cr  -f  Cp  <  0  and  Cp  -f-  2Ch  <  0,  both  polarizations  are  evidently  modu¬ 
lationally  stable,  contrary  to  the  claim  of  Chian  and  Kennel  [1].  As  noted  above,  their 
solution  violates  the  z-component  of  the  force  equation  at  leading  order  and  thus  ignores 
ponderomotive  effects,  which  (as  we  saw)  are  siabilizing  for  both  polarizations.  For  linear 
polarization,  the  additional  stabilizing  effect  of  mode-mode  interactions  (harmonic  effect) 
was  apparently  also  neglected. 

If  the  electrons  and  positrons  have  finite  (but  equal)  temperatures,  we  find 


Cp  +  2  Cr  = 


2  2 
e  Up 

m2c* 


(29c) 
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Cp  +  3  Cr  +  Ch 


m2c4 


e2u>p 

m2c4 


1  ("l-c2±k>\  _  (  ul  +  cl/cW  \ 

2  yu2  —  c±k2  J  \^c2k2  —  c2±[ c2ui2  J 

1  ( (1  +  c -Jc2)u>l  -  c2±/c2lv2\ 

2  y(l  -  c2±/c2)u2  +  c|/c2u>2  J 

(  ul  +  c2±/c2v2  \ 

^2(l-^-)-u ,j). 


(29/) 


where  c±  is  the  speed  of  sound  of  both  components.  Due  to  the  singularity  in  the  pon- 
deromotive  coefficient,  modulational  instability  in  a  electron-positron  plasma  is  possible, 
but  only  in  a  narrow  range  just  :  bove  the  plasma  frequency. 


We  now  briefly  discuss  astrophysical  implications:  Recent  observational  data  [14-16] 
strongly  support  the  hypothesis  that  pulsar  micropulses  are  a  ” temporal”  phenomenon  and 
can  be  interpreted  within  the  amplitude-modulated  noise  model  [17].  As  discussed  in  [2-5], 
nonlinear  modulational  instability  would  provide  a  natural  mechanism  for  amplitude  mod¬ 
ulations.  Now,  according  to  pulsar  models  [18-19],  the  pulsar  magnetosphere  is  composed 
of  secondary  electrons  and  positrons.  However,  this  paper  shows  that,  in  the  absense  of 
an  ambient  magnetic  field,  a  cold  electron-positron  plasma  does  not  exhibit  modulational 
instability  for  linear  or  circular  polarization  and  that  finite  temperature  affects  this  result 
only  just  above  the  plasma  frequency.  On  the  other  hand,  we  have  found  that  a  cold 
electron  positive-ion  plasma  is  mod  Nationally  unstable  for  either  polarization.  So  if  pulsar 
micropulses  really  are  caused  by  modulational  instability,  then  at  least  one  new  feature 
must  enter  the  problem. 

We  would  like  to  mention  three  possible  candidates  for  this  new  feature  and  of  interest 
for  pulsar  micropulses:  First,  from  our  above  result  on  the  modulational  instability  of  an 
electron-positive  ion  plasma,  one  would  expect  that  the  inclusion  of  some  positive  ions  coui  \ 
drive  the  electron-positron  plasma  modulationally  unstable.  In  a  three-component  neutra1 
plasma  consisting  of  electrons,  positrons,  and  positive  ions,  we  expect  that  a  critical  ratic 
(positive-ion  density)/(positron  density)  will  exist  -  depending  perhaps  on  parameters  such 
as  the  frequency  -  since  the  limiting  cases  (electron-positron,  electron-ion)  give  opposing 
results. 


The  second  physical  effect  would  be  the  presence  of  ambient  magnetic  fields,  which 
could  strongly  affect  the  modulational  instability  of  any  plasma,  in  particular  an  electron- 
positron  plasma.  Careful  calculations  are  required  in  order  to  determine  how  such  a  mag¬ 
netic  field  would  affect  the  nonlinear  coefficients,  keeping  in  mind  that  new  parametric 
interactions  can  occur  due  to  the  presence  of  the  ambient  magnetic  field.  (Such  a  cal¬ 
culation  was  proposed  in  [5],  but  unfortunately  the  treatment  contains  errors,  as  can  be 
confirmed  for  example  by  comparison  of  Eq.  (7)  of  [5]  with  Eq  (44)  of  [20].) 

Third,  in  pulsar  magnetospheres,  the  actual  estimated  value  of  the  strength  parameter 
of  Eq.  (2)  may  be  of  order  unity  or  larger.  Under  these  circumstances,  the  present  theory 
cannot  be  applied  directly,  since  higher-order  corrections  to  the  nonlinear  coefficients  found 
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here  could  be  important,  and  at  the  moment,  no  a  priori  method  for  estimating  their 
size  is  available.  A  fully  relativistic  scheme  allowing  arbitrary  (slow)  modulations  of  a 
zeroth-order  solution  with  large  (relativistic)  amplitudes  would  be  optimal.  However,  such 
a  scheme  will  contain  all  the  subtleties  of  the  present  case  together  with  possibly  new 
ones.  (For  example,  slowly  varying  “DC”  ponderomotive  terms  in  the  longitudinal  force 
equations  will  enter  two  orders  in  e  sooner  than  in  the  present  calculation.)  Now,  as  we 
saw  above,  the  Karpman-Krushkal  [6]  application  of  Whitham  theory  [7]  as  applied  in  [1] 
and  [3-5]  unfortunately  did  not  give  satisfactory  results  in  the  present  problem.  But  we  do 
note  that  if  the  Whitham  method  could  be  extended  to  include  parametric  interactions, 
then  it  might  w'ell  be  possible  to  obtain  estimates  of  nonlinear  coefficients  even  in  the  fully 
relativistic  case. 
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Abstract  \ 

This  paper  extends  our  previous  results  [1]  on  the  nonlinear  modulational  stability 
properties  of  plasma  electromagnetic  pulses,  to  include  the  presence  of  an  ambient  mag¬ 
netic  field  Bo  parallel  to  the  direction  of  propagation.-  As  before,  the  pulse  is  assumed 
to  be  strong  enough  to  accelerate  particles  to  weakly,  but  not  fully,  relativistic  velocities. 
The  positive  component  may  consist  of  either  positrons  or  singularly  charged  ions,  and  no 
specific  assumptions  or  approximations  are  made  concerning  the  mass  ratio  of  the  com¬ 
ponents.  (Previous  work  assumed  a  positron-electron  plasma.)  The  plasma  is  assumed 
fully  ionized.  The  effects  of  a  finite  temperature  are  included  for  generality.  Using  singular 
perturbations,  we  derive' approximate  solutions, which  describe  the  evolution  of  a  circularly 
polarized  pulse.  The  envelope  is  shown  to  evolve  over  long  time  scales  according  to  the 
cubic  nonlinear  Schroedinger  (NLS)  equation.  Relativistic  corrections  and  ponderomotive 
forces  both  contribute  terms  cubic  in  the  amplitude.  (In  contrast  to  the  case  studied  in 
[1],  harmonic  effects  vanish  identically  here  because  of  circular  polarization.)  A  positron- 
electron  plasma  without  magnetic  fields  was  shown  in  our  previous  paper  to  be  modula- 
tionally  stable,  except  in  the  case  of  finite  temperature  where  modulational  instability  is 
possible  near  the  plasma  frequency  uv.  Here,  it  is  shown  that,  even  for  a  cold  plasma,  the 
presence  of  an  ambient  magnetic  field  makes  a  decisive  difference:  Modulational  instability 
can  arise  within  a  broad  range  of  frequencies  and  values  of  JE?o,  in  particular  for  a  pure 
positron-electron  plasma.  For  given  Bo  and  polarization,  we  demonstrate  the  existence 
of  critical  frequencies  for  the  onset  of  modulational  instability.  This  result  has  important 
consequences  for  observations  of  pulsar  micropulses  and  possible  technological  applications. 
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The  possible  influence  of  ambient  magnetic  fields  on  the  modulationa!  (in)stability  of 
electromagnetic  pulses  in  weakly  relativistic  plasmas  has  taken  on  new  importance  with 
the  recent  finding  [1]  that  a  cold  positron-electron  plasma  without  ambient  magnetic  fields 
is  modulationally  stable  for  both  linear  and  circular  polarization,  contrary  to  previous 
findings  [2-5].  However,  there  is  strong  recent  observational  evidence  [6-9]  that  pulsar  mi¬ 
crostructure  may  result  from  modulations  of  a  coherent  pulse.  If  nonlinear  modulational 
instability  is  the  cause  of  micropulses,  then  some  additional  element  must  enter  the  pic¬ 
ture:  e.g.,  finite  temperature,  ambient  magnetic  fields,  possible  presence  of  ions,  or  fully 
relativistic  amplitudes.  As  shown  in  [1],  finite  temperature  can  indeed  induce  modulational 
instability  for  both  circular  and  linear  polarization,  but  the  effect  is  confined  to  frequen¬ 
cies  near  the  plasma  frequency.  In  this  paper,  we  show  that  in  the  presence  of  an  ambient 
magnetic  field,  a  circularly  polarized  pulse  propagating  parallel  to  the  magnetic  field  in  a 
cold  plasma  can  go  modulationally  unstable  under  a  range  of  conditions. 

A  previous  attempt  [5]  to  determine  conditions  for  soliton  solutions  in  a  strongly 
magnetized  plasma  began  from  a  linear  solution  which  unfortunately  contained  errors,  as 
pointed  out  in  [1].  These  errors  were  partially  corrected  in  [10].  However,  the  results 
apply  only  to  the  special  case  of  a  single-soliton  solution.  Correct  formulas  for  the  linear 
theory  can  be  found  in  [11-12].  As  applied  in  [2-5],  the  formulas  given  by  Karpman  and 
Krushkal  [13]  (based  on  Whitham  theory  [14])  resulted  in  expressions  which  do  not  satisfy 
the  Maxwell  equations  and  the  Lorentz  force  law,  and  the  nonlinear  coefficients  of  the  NLS 
equation  were  as  a  result  incorrect. 

For  these  reasons,  we  have  applied  the  singular  perturbs*' on  method  (“two-tirrung”) 
used  in  [1]  for  the  purpose  of  accurately  evaluating  the  nonlinear  coefficnets.  (For  economy, 
we  refer  the  reader  to  the  explanations  given  there.)  Since  the  physical  size  of  a  term  is 
explicitly  taken  into  account  in  assigning  orders  of  smallness,  the  method  contains  internal 
checks  which  do  not  allow  one  to  ’’lose”  terms.  We  have  employed  the  computer  algebra 
system  MUMATH  to  verify  that  all  equations  are  indeed  satisfied  at  the  orders  claimed. 
(The  two-timing  techniques  of  [1]  were  automated  by  rewriting  the  usual  differentiation 
routines.)  References  to  equation  numbers  from  [1]  in  the  following  will  be  preceeded  by 
Roman  numeral  I. 

Following  [1],  our  model  comprises  a  singly  ionized,  positively  charged  species  of  mass 
M  and  a  negatively  charged  species  of  mass  m,  hereafter  refered  to  as  the  ions  and  electrons, 
respectively.  (However,  since  no  approximation  is  made  at  this  stage  with  respect  to 
m/M,  all  formulas  will  be  equally  valid  for  an  electron-positron  plasma  upon  substituting 
M  =  m.) 

The  fluid  is  characterized  as  in  [1]  by  the  electron  density  n,  the  electron  velocity  v, 
the  ion  density  N,  the  ion  velocity  V  the  vector  potential  A  and  the  electrostatic  potential 
The  function  n,  v,  N ,  F,  and  <f>  as  well  as  the  perturbed  part  of  the  vector  potential 
A  are  assumed  to  depend  on  t  and  z,  but  not  on  x  and  y.  The  unperturbed  part  of  A 
describes  a  uniform  magnetic  field  in  the  z-direction  and  thus  of  course  depends  linearly 
on  x  and  y.) 
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In  addition  to  the  plasma  frequency,  given  by 


= 47re2no  u + 

(1) 

we  define  the  electron-cyclotron  frequency 

u>-  =  e|B0|/(mc) 

(2a) 

and  the  ion-cyclotron  frequency 

w+  =  e|B0|/(Mc) 

(26) 

where  c  is  the  speed  of  light.  (The  two  cyclotron  frequencies  obviously  coincide  in  the  case 
of  an  electron-positron  plasma.) 

As  in  [1],  we  suppose  that  the  transverse  electric  field  E±  of  the  electromagnetic  pulse 
satisfies 

e\Ex\/(mc)  =  0(e),  (3) 

where  e  is  a  dimensionless  small  parameter  and  e  <<  1.  As  a  consequence,  the  relativistic 
factors  7  and  T  (defined  below)  differ  from  unity  by  0(e2). 

It  is  convenient  to  express  the  continuity  and  Euler  equations  for  the  electrons  in  the 

form 

dtn  +  V  •  (nv)  —  0  (4 a) 

d(myv  —  -A)  —  eV<f>  +  me,2  — Vn  +  -E^=1UfVAf  =  0  (46) 

c  n  c 

where 

d  =  dt  v  •  V  (4c) 

7  =  (1  -  n2/c2)-1/2  (4 d) 

and  c,  is  the  electron  thermal  speed.  For  the  ions,  the  respective  equations  are 

dtN  +  V  ■  (NV)  =  0  (5a) 

D(MTV  +  -A)  +  eV<f>  +  -  -X3t=1VtVAt  =  0  (56) 

c  Ac 

where 

D  =  dt  +  V  ■  V  (5c) 

r  =  (1  -  V2/c2)~1/2  (5 d) 

and  C,  is  the  ion  thermal  speed.  Of  course,  Maxwell’s  equations 

_  _  —  1 

VxB  =  - (NV  -nv)+-dtE 

c  c 

VE  =  4ne(N  -  n) 
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(ба) 

(бб) 


must  hold,  where  we  take 


5  =  VxI 

(6c) 

E  =  -V<?$  -  -dtA 
c 

(6  d) 

and  where  we  impose  the  radiation  gauge 

o 

II 

»o 

(6e) 

An  exact  static  solution  of  the  system  (4-6)  is  given  by 

3 

II 

II 

3 

o 

(7a) 

v  =  0 

(7b) 

V  =  0 

(7c) 

o 

il 

(7  d) 

A  =  A0  =  (— yB0/2,  iJ3q/2,0) 

(7c) 

which  contains  a  constant  magnetic  field  in  the  z  direction.  Quantities  with  subscript  0 

are  strictly  constant. 

Following  [1],  we  assume  asymptotic  expansions  of  the  form 

n  =  no  +  enj  +  e2n2  +  e3n3  +  . . . 

(8a) 

A  =  A0  +  eA\  +  e2  A2  +  e3A3  +  . . . 

(8b) 

N  =  n0  +  eN1+  e2JV2  +  e3  N3  +  ... 

(8c) 

v  =  ev[  4-  e2V2  +  e3vl  -f  . . . 

(8d) 

V=  eV1+e2V2  +  e3V3  +  ... 

(Be) 

<f>  =  e<j>  i  A  e2  4>2  e3  <f>3  +  . . . 

(8/) 

where  Aq  is  given  by  (7e). 

It  is  instructive  to  review  briefly  properties  of  the  solutions  of  the  linearized  system 
obtained  by  simply  truncating  Lqs.  (4-6)  at  0(e).  (The  linearized  theory  of  circularly- 
polarized  electromagnetic-wave  propagation  parallel  to  an  ambient  magnetic  field  can  be 
found  in  [12].)  A  typical  Fourier  component  of  the  vector  potential  perturbation  behaves 
in  this  linearized  theory  like  a( u>)e'(k*~wi\  where 

a  =  a-(x  -f  ty).  (9) 

The  dispersion  relation  for  right-hand  circular  polarization  is 

u)2  =  47re2no/(mr')  +  47re2n0/(A//r)  +  k2c2  (10) 
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where 


v  =  1  —  eBo/irncui)  =  1  —  u-/u>  (Ha) 

fi  =  1  +  eB0/(Mcu)  =  1  +  u>+/u>  (116) 

The  dispersion  relation  for  left-hand  circular  polarization  can  be  obtained  from  (11)  by 
substituting  for  v  and  /j.  expressions  in  which  Bo  is  replaced  by  —Bo. 

If  the  masses  are  unequal,  then  four  electromagnetic  branches  can  be  distinguished. 
If  the  masses  are  equal  (electron-positron  plasma),  then  two  branches  are  present.  The 
requirements  k2  >  0  and  u>2  >  0  will  in  general  place  restrictions  on  the  allowable  frequen¬ 
cies  of  a  propagating  electromagnetic  wave.  Unlike  the  unmagnetized  case,  propagation 
can  occur  below  the  plasma  frequency.  In  the  linear  theory  it  is  of  course  also  possible 
to  consider  the  evolution  of  initial  data  corresponding  to  “linear”  polarization;  evidently 
Faraday  rotation  of  the  polarization  will  occur. 

In  addition  to  the  electromagnetic  modes,  there  are  also  longitudinal  Langmuir  oscil¬ 
lations  at  the  plasma  frequency.  In  the  linearized  theory,  these  are  of  course  decoupled 
from  electromagnetic  waves.  However,  in  the  modulated  nonlinear  theory  we  shall  see  that 
the  ponderomotive  force  induces  coupling  between  transverse  and  longitudinal  modes. 

Following  the  procedure  of  [1],  we  now  assume  that  the  first-order  perturbation  of  the 
vector  potential  and  electron  density  take  the  form 

Ai  =  &«**—*>  +  c.c.  (12) 

a  -  a  ■  (x  +  iy)  (13) 

ni=0,  (14) 

where  ”c.c.”  indicates  the  complex  conjugate,  and  where  the  wave  number  k  and  the  fre¬ 
quency  u  correspond  to  nonsingular  solutions  of  the  dispersion  relation  (10).  In  particular 
we  assume  that  u>  is  quite  distinct  from  the  electron-cyclotron  frequency  and  also  the  ion- 
cyclotron  frequency.  (In  fact,  we  exclude  a  small  range  of  0(e)  around  these  frequencies). 
The  function  a  is  not  strictly  constant,  but  depends  on  the  slow  variables 

Z  =  ez,  T  =  et,  t  =  e2t  (15) 

This  choice  of  the  form  of  the  first-order  perturbation  places  restrictions  on  the  class  of 
initial  conditions  described  by  our  treatment:  The  initial  conditions  correspond  to  a  right- 
hand  circularly  polarized  pulse  propagating  in  the  positive  z  direction.  (Results  for  the 
opposite  polarization  can  be  recovered  by  replacing  Bo  with  —Bo  at  the  end.) 

As  in  [1],  in  first  order,  the  longitudinal  components  are  identically  zero  or  can  be 
chosen  to  vanish  without  loss  of  generality.  The  transverse  components  can  be  treated  as 
a  coupled  linear  homogeneous  system  of  the  form 

L(rP  i)  =  0  .  (16) 
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The  solution  can  be  expressed  in  the  form 


(17a) 

(176) 

(17c) 


4>\  =  0  =  ni  =  Ni  =  Ait  =  =  Viz 

i?i  =  e^iz-w^ea/(mct/)  +  c.c. 

V!  =  -e,(il'ul)ea/(Mc/i)  +  c.c. 

The  notation  means  the  the  0(e)  part  of  v,  and  so  forth.  The  dispersion  relation  (10)  of 
the  (truncated)  linearized  theory  must  also  hold.  The  first-order  part  of  (6b)  guarantees 
quasineutrality  at  first  order.  The  2-component  of  (6a)  is  then  satisfied  at  first  order 
automatically. 

We  now  refer  the  reader  to  Eqs.  (1. 15)-(1. 1 8)  and  discussion.  As  before,  the  second- 
order  parts  of  the  transverse  equations  -  together  with  the  condition  for  no  secular  terms 
-  restrict  the  slow  dependence  of  a.  One  finds  explicitly  that  the  amplitude  is  transported 
according  to 


(&r  4-  vgdz )  a  =  0 

(18a) 

at  the  group  velocity 

vg  =  dui/dk 

(186) 

where 

II 

P 

/  [l/(mV)_  1/(A,V)]  ) 

(IS  c) 

where  v  and  p  were  defined  in  Eq.  (11)-  Note  that  the  product  of  group  and  phase  velocity 
is  not  c2.  We  now  examine  the  next  leading  order  in  the  remaining  equations.  These  consist 
of  third  order  in  the  longitudinal  parts  of  the  Euler  equations  (4b, 5b),  second  order  in  the 
continuity  equations  (4a, 5a),  second  order  in  Poisson’s  equation  (6b),  and  second  order 
in  the  transverse  Euler  equations  (4b, 5b).  (The  longitudinal  part  of  (6a)  is  then  satisfied 
automatically.) 

We  recall  that  in  [1],  the  second-order  solution  contained  three  types  of  terms  classified 
according  to  their  phase  dependence  as  “DC,”  “fundamental,”  and  “harmonic.”  All  three 
types  of  terms  contributed  effects  of  comparable  magnitude,  in  the  general  case.  However, 
harmonic  terms  vanished  in  the  special  case  of  circular  polarization.  The  same  is  true  here, 
even  though  a  magnetic  field  is  present.  On  the  other  hand,  ponderomotive  forces  arise  in 
the  longitudinal  equations  and  induce  DC  corrections  as  before. 

The  second-order  equations  form  a  linear  inhomogeneous  system,  and  thus  given  any 
particular  second-order  solution  one  can  find  another  particular  solution  (corresponding 
to  different  initial  data,  of  course)  by  the  addition  of  an  arbitrary  homogeneous  solutions 
of  the  form  (16). 

Ai  =  0,  (19) 

we  obtain  for  the  DC  part 


ri2  =  Q.Q*  TIqO 


(20a) 
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where 


AT2  =  ri2 
V2  z  =  vgaa*a 
V^z  =  vgaa*<7 
<j>2  —  a.a*Tj 


e2(u>2  4  c2k 2  —  2kvgu> ) 
c2us2mMvg2  F 


(206) 

(20c) 

(orW) 

(20e) 


(21) 


£  = 


(Af  +  m)c2vg2F 


— (4-c,2) 

mi/ 


1  + 


kvgu>c 

u2u 


m 


Mu 


(V-c.2) 


1  - 


kVgQc 
uj2  H 


(22) 


me,2  4  MC,2 
vg7(j\I  4  m) 


(23) 


Note  that  F  contains  the  corrections  due  to  finite  temperature.  Also  note  that  longi¬ 
tudinal  perturbations  are  forced  by  the  modulations  in  a.  As  in  [1],  charge  quasi-neutrality 
at  second  order  still  holds.  However,  in  contrast  to  [1]  the  second-order  transverse  velocity 
corrections  are  nonvanishing.  These  are 


v2x  =  -ie,(k*  w<)e2So^rd/(mzu.-zv2)  4  c.c.  (24a) 

V^jl  =  -i :.eilkt'u,,)e2BodTa/(M2u}2n2)  4  c.c.  (246) 


We  now  proceed  to  third  order.  Choosing  the  initial  conditions  as  above  to  exclude 
homogeneous  Langmuir  oscillations  and  again  absorbing  A3  into  the  first  order,  we  solve 
the  transverse  parts  of  the  Euler  equations  (4b, 5b)  to  obtain  the  particular  solution 


V3x  =  e,(kz  “f)(i  4  ty)s3  4  c.c. 

V3_l  =  °(f  4ty)534c.c. 

where  s3  is  given  by 

s3  =  —  ieu_(9ra/(mcw’2i/2)  -  eu >..dj-a/(mcvs3 1/3) 

4  ekvgu>-a2a‘o/(ma. ?v2)  -  2 e3aza*/(m3i/4) 


(25a) 

(256) 


(26) 


An  expression  for  53  can  be  obtained  from  _s3  by  replacing  m  with  M,  v  with  fi,  u'_  with 
w+ ,  and  e  with  -e.  The  last  term  in  (26)  is  the  relativistic  correction  and  except  for  the 
corresponding  term  in  S3,  it  is  the  only  relativistic  correction.  The  other  nonlinear  term 
in  (26)  is  a  ponderomotive-like  term  arising  from  the  second-order  shift  in  the  longitudinal 
velocity,  vzx. 
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The  longitudinal  equations  at  the  next  leading  order  give  information  about  the  third- 
order  density  and  potential  corrections.  However,  we  already  have  all  the  information 
necessary  to  evaluate  the  secular  condition  implied  by  the  transverse  components  of  the 
Maxwell  equation  (6a).  The  linear  terms  are  consistent  with  what  would  have  been  ob¬ 
tained  by  the  Karpman-Krushkal  [13]  approach.  In  particular,  the  coefficient  of  dTa  can 
be  shown  to  be  equal  to  2 kc2/vg)  The  dispersion  terms  (second  derivatives  in  T  and/or  Z) 
can  be  collected  together  by  means  of  (18a)  into  a  coefficient  times  d\a.  This  coefficient 
turns  out  to  be  kc^Ukk/vg,  where  aytt  is  called  the  dispersion  and  is  the  second  derivative 
of  u  with  respect  to  k.  To  be  sure,  considerable  algebra  is  required  in  order  to  obtain  the 
final  result.  However,  there  are  some  ’’miraculous  cancellations”  and  things  do  simplify 
toward  the  end.  After  division  by  some  common  factors,  we  finally  arrive  at  an  equation 
of  the  form 

2iudra  -f-  uju>kk9^Q  ~h  a2a  (C p  -h  2 C  r)  —  0  (2  i) 

where  the  nonlinear  coefficient  Cr  contains  the  effects  of  the  relativistic  corrections,  while 
the  nonlinear  coefficient  Cp  contains  the  ponderomotive  effects.  The  latter  arises  from 
two  sources.  Firsi  there  is  the  shift  in  the  second-order  longintudinal  velocity  mentioned 
above.  This  term  did  not  occur  in  the  zero-magnetic-field  case  [1],  The  other  source  is 
the  second-order  shift  in  the  density  which  did  occur  in  the  zero-magnetic-field  case  [1], 
Finally,  note  that  no  terms  of  the  type  refered  to  in  [1]  as  ’’harmonic”  occur  here  (such 
terms  would  be  expected  for  generic  polarization). 

The  sum  of  the  two  ponderomotive  terms  is  given  simply  by 


Cp  -  - 


(2  So) 


The  relativistic  coefficient  is  giver,  explicitly  by 

,2 


Cr  = 


uJVQM2m2u.'l  r  ,  ,  .  ,  .  , 

— - MV(mV)  +  l/(A/V)] 


c4A/m  /  2c2k(M  -I-  m) 


;2S  b 


We  note  that  the  coefficients  Cp  and  Cr  reduce  to  the  corresponding  formulas  (I.25b-c)  if 
B o  vanishes. 

The  condition  for  modulational  instability  is 


vukk{Cp  +  2 Cr)  <  0 


(29) 


It  is  remarkable  that  Cr  can  in  fact  be  reduced  to  a  functional  of  the  linear  dispersion 
relation  alone.  One  can  show  that  (28b)  is  equivalent  to 


Cr 


u>  va 


2c2£u>_u>+ 


W(c2*J)-l  +  iay(rV) 
6  2 


(30 


where  k  is  regarded  here  as  an  implicit  function  to  be  obtained  by  solving  the  dispersion 
relation  for  k  in  terms  of  us.  The  practk.d  advantage  of  (30)  (together  with  (2Sal)  is  that  it 
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reduces  the  problem  of  evaluating  the  NLS  coefficents  to  formulas  involving  only  the  linear 
dispersion  relation.  Given  c2k2(u> ),  one  may  now  evaluate  all  of  the  NLS  coefficients.  Thus 
the  problem  of  calculating  the  nonlinear  coefficents  is  simply  a  problem  in  evaluating  vg 
and  Slj{c2k2).  These  expressions  could  be  useful  in  obtaining  clues  to  the  corresponding 
result  for  arbitrary  polarization. 


We  examine  the  electron-positron  case  in  detail.  In  this  case,  the  dispersion  relation 
(10)  reduces  to 

.2  T 


c2k2  =  u,2 


(jJ 

1 _ 


w2  -  n2 


(31) 


where  ft  is  now  the  common  cyclotron  frequency.  The  group  velocity  is  found  to  be 


u  )  ,  . 

'  1  +  TP-TF? 


and  the  nonlinear  coefficients  are 

.2 


c*=\74 


U2p^Vg 


c4m2  J  2/cc2(u.'2  -  O2)2 


c4m2  )  Jfc3c4(w2  -ft2)2F 


1  + 


8u2Q.2 

1/  +  (u>*  -  n2)2  J 

2ft2 


a>2ft2 


U>2-ft2  (u>2  —  ft2  )2 


while  the  dispersion  is 


Vkk  = 


v*u>2u>2 

9  P 


c4fc3(w2  —  ft2) 


1  + 


5ft2  ft2(4ft2  —  3u;2)  3u;2ft4 


+ 


(u>2-ft2)  (u>2-ft2)2  (u>2-ft2)3 


(32) 


(зза) 

(ззб) 


(34) 


In  the  limit  of  large  ft2/w2,  which  is  certainly  a  good  approximation  for  the  pulsar  envi¬ 
ronment,  one  obtains  from  the  above,  for  kjus  >  0, 


C*  = 


Cp  =  - 


=c(l-^+---) 


Vkk  =  - 


2ft2 
3o >u2c2 


ft4 


+ 


c4m2 

„2 


c4m2 


u>p2u)ivg 
2fcc2ft4 

U>p2U j*Vg 

kc2U4F 

-2  +  C? 


1  + 


10u>2 


+ 


1  + 


ft2 

6u>2  +  u>„2 


ft2 


F  =  1  - 


2c2 


+ 


0 


+  •••) 


(35a) 

(356) 
(35  c) 

(35c) 


(35c) 
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Note  that  (35b)  gives  u>kkw  <  0.  In  the  limit  of  Bq  — »  0,  we  have  the  opposite 
result  of  u>k *<*>  >  0.  Recall  that  in  I,  it  was  positive  nonlinear  coefficients  which  led  to 
modulational  instability;  here  it  will  be  negative  nonlinear  coefficients  which  will  give  rise 
to  modulational  instability. 

As  long  as  the  sound  speeds  are  much  smaller  than  the  group  velocity  of  the  electro¬ 
magnetic  wave,  (which  is  usually  the  case  for  nonrelativistic  gas  temperatures)  F  is  close 
to  unity.  (However,  the  A  and  FMS  branches  [12]  of  the  dispersion  relation  have  van¬ 
ishing  group  velocity  near  the  cyclotron  frequencies,  and  therefore  for  finite  temperature 
sufficiently  near  these  frequencies  F  may  become  significantly  different  from  unity.)  Using 
Eqs.  (35),  we  find 


Cp  4-  2  Cr 


\c2m2  / 


u>p2u4 

2ft4 


4  u>2  —  u>l 

_ f 

ft2 


i  ±  £ ,2 

2c2 


+ 


(36) 


Now  we  can  have  modulational  instability,  see  (29),  only  if 


2  1  2  r,2Cl+Cs 

u>  <  -u  -(-  ft2  — - - 

4  r  2c2 


(37) 


For  large  ft2,  we  see  that  there  is  a  very  wide  frequency  range  for  even  moderate 

temperatures.  And  even  for  a  cold  plasma,  modulational  instability  occurs  also  if  u2  < 
i,  ,  2 
4UP  ■ 

According  to  pulsar  models  [15,16],  the  pulsar  magnetosphere  is  composed  of  sec¬ 
ondary  electrons  and  positrons,  and  strong  magnetic  fields  (1012  Gauss)  are  present,  cor¬ 
responding  to  an  electron  plasma  frequency  of  the  order  1019  Hz,  as  compared  to  a  plasma 
frequency  of  a  few  Megahertz.  If  we  apply  the  statistical  model  of  Rickett  [17],  as  described 
in  [2],  a  modulationally  unstable  pulse  could  be  responsible  for  the  micropulse  structure 
observed  in  [18]. 

Unfortunately,  one  cannot  yet  make  any  definitive  predictions  for  the  problem  of 
micropulses.  First,  one  would  expect  that  the  propagation  angle  will  in  general  not  be 
parallel  to  the  magnetic  field.  Second,  the  polarization  may  not  necessarily  be  circular. 
Third,  the  pulses  may  well  be  fully  relativistic.  Still,  based  on  the  present  calculation, 
one  would  expect  that  a  range  of  frequencies  for  modulational  instability  should  exist 
even  when  the  propagation  angle  is  not  parallel  to  the  magnetic  field,  and  even  when  the 
polarization  is  elliptic.  (Perhaps  the  observed  mixed  polarizations  result  from  averaging 
due  for  example  to  insufficiently  fine  time  resolution.)  Consideration  of  fully  relativistic 
effects  will  have  to  await  the  results  of  further  work. 


Our  paper  predicts  a  definite  frequency  range  for  instability;  for  the  case  of  a  positron- 
electron  plasma,  the  frequency  range  is  given  directly  in  Eq.  (37).  For  the  more  general 
case,  the  frequency  range  can  easily  be  obtained  from  Eqs.  (28)  and  (29),  together  with 
the  dispersion  relation  (10). 
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A  further  result  of  our  calculation  is  a  prescription  for  obtaining  the  nonlinear  coeffi¬ 
cients  directly  from  the  linear  dispersion  relation,  as  in  Eqs.  (30)  and  (28a)  above.  Note 
that  these  two  equations,  although  obtained  for  a  two  component  plasma,  have  an  obvi¬ 
ous  generalization  for  a  plasma  containing  any  number  of  components.  Thus,  given  an 
N-component  neutral  plasma  and  a  linear  dispersion  relation,  one  can  directly  obtain  the 
coefficients  of  nonlinearity  in  the  case  of  circular  polarization.  This  is  in  the  spirit  of  the 
Karpman-Krushkal  [13]  method,  which  Refs.  [2-5]  and  [10]  attempted  to  apply. 
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